*the following is the fraction of persons with occupations that earn lower income than the 'vulnerability' level as defined in the paper
replace vulnerable=0.6088632 if ct=="bdf"
replace vulnerable=0.7167798 if ct=="bek"
replace vulnerable=0.5765517 if ct=="buk"
replace vulnerable=0.660061 if ct=="cam"
replace vulnerable=0.5105602 if ct=="chs"
replace vulnerable=0.3596974 if ct=="con"
replace vulnerable=0.3831878 if ct=="cul"
replace vulnerable=0.4966975 if ct=="dby"
replace vulnerable=0.4060673 if ct=="dev"
replace vulnerable=0.6514936 if ct=="dor"
replace vulnerable=0.3222106 if ct=="dur"
replace vulnerable=0.5902228 if ct=="ess"
replace vulnerable=0.5209353 if ct=="gls"
replace vulnerable=0.589928 if ct=="ham"
replace vulnerable=0.6791745 if ct=="her"
replace vulnerable=0.6860215 if ct=="hrt"
replace vulnerable=0.7492877 if ct=="hun"
replace vulnerable=0.5342939 if ct=="ken"
replace vulnerable=0.4559278 if ct=="lan"
replace vulnerable=0.5450517 if ct=="lei"
replace vulnerable=0.670229 if ct=="lin"
replace vulnerable=0.3699086 if ct=="mdx"
replace vulnerable=0.4246753 if ct=="nbl"
replace vulnerable=0.6004073 if ct=="nfk"
replace vulnerable=0.5594901 if ct=="nth"
replace vulnerable=0.4821853 if ct=="ntt"
replace vulnerable=0.7155655 if ct=="oxf"
replace vulnerable=0.5892857 if ct=="rut"
replace vulnerable=0.56089 if ct=="sal"
replace vulnerable=0.5242395 if ct=="som"
replace vulnerable=0.3657289 if ct=="sta"
replace vulnerable=0.5709323 if ct=="suf"
replace vulnerable=0.5515643 if ct=="sur"
replace vulnerable=0.6170678 if ct=="sus"
replace vulnerable=0.4659166 if ct=="war"
replace vulnerable=0.6518595 if ct=="wil"
replace vulnerable=0.4692737 if ct=="wln"
replace vulnerable=0.4703344 if ct=="wor"
replace vulnerable=0.4689523 if ct=="yks"
